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We calculate massive 5-propagator 2-loop integrals for operator matrix elements in the light-cone expansion, using 
Mellin-Barnes techniques and representations through generalized hypergeometric functions. 
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1. Introduction 

The use of Mellin-Barnes integrals became a 
widespread technique for calculating Feynman di- 
agrams throughout the last years in partic- 
ular to calculate double and triple box-diagrams. 
In Ref. [12], it was possible to expand the scalar 
two-loop two-point function in all orders in the 
dimensional regularization parameter e, using ad- 
ditionally the gluing operation of Feynman dia- 
grams, defined by Kreimer [0]. In this paper, 
we will apply this technique to a more complex 
problem, namely the calculation of massive five- 
propagator 2-loop Feynman Diagrams with oper- 
ator insertions, stemming from light-cone expan- 
sion, which are needed for the calculation of the 
heavy flavor coefficient functions in deep-inelastic 
scattering. These Wilson coefficients have been 
calculated before up to next-to- leading order 
Fully analytic results could only be obtained in 
the limit Q 2 3> m 2 using mass factorization 

The heavy flavor Wilson coefficients are ob- 
tained as a convolution of the massless Wil- 
(Q 2 / M 2 ? z ) an d the massive 



son coefficients C^ L 



operator matrix elements (OME) A\ J (/n 2 /m ), 
which shall be calculated to order 0(a 2 s ). A cal- 
culation of these operator matrix elements was 
performed in [|3] using integration-by-part tech- 
niques [ Ej, which leads to a large number of 
terms. Expressing the result in terms of Nielsen 
integrals leads to lengthy expressions for most in- 
dividual diagrams and the complete result. 

This paper offers a different approach, using 
Mellin-Barnes integrals, which will enable us to 
obtain more simple results, with the positive side- 



effect, that we avoid the creation of many of the 
above-mentioned canceling intermediate terms. 
More precisely, we will directly obtain analytic 
results in Mellin space, for a general value of the 
Mellin variable N, using multiple harmonic sums 
[ Ej , which allows to compactify the result ob- 
tained in x-space [EJ significantly. Before we will 
do this, we briefly review the method used in [Ej 
and starting from this proceed to the more com- 
plex cases at hand. Finally, we will present results 
obtained in this way. 

2. The Method 

For the process we want to calculate, we will en- 
counter diagrams with operator insertions located 
cither at a line or a vertex of the diagram. In 
Figure 1 the seven five-propagator integrals con- 
tributing arc shown. The powers vi of propaga- 
tors are one, except for diagram I, where also the 
power v\ = 2 occurs for the line with the opera- 
tor insertion. These diagrams will be calculated 
in a similar way to the calculation performed for 
the massless two-loop two-point function, cf. 
12 IB]- We start here with the massive two- loop 
two-point function, with four massive propaga- 
tors (thick lines). The two-loop two-point dia- 
gram can be split into the one-loop two-point and 
the one-loop three-point function, using the glu- 
ing operation of graphs in D = 4 — 2e dimensions 
[0: 
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Figure 1. 
sertions; 



The seven diagrams with operator in- 
-®- = {A. P ) N -\ 



gram is re-obtained by inserting the three-point 
function into the two-point function. The same 
as for the graphs is done on the level of inte- 
grals, "splitting" the two-loop integral into its 
one-loop two-point and three-point components. 
The insertion of the three-point function into the 
two-point function is reflected by the fact that 
the result of the three-point function alters the 
result of the two-point function by shifting the 
arguments of its propagator-exponents. We start 
with the three-point function: 



J [k 2 + m 2 }^ pa" - p) 2 + m 2 ]"3 [(fc 2 - fc^ 2 ]^ 

In a first step, we apply Feynman parameteriza- 
tion to the propagators, turning the product of 



three propagators into a sum of them: 



1 



■sin 



dx\...dx n 5 



(5> 



ILx? 



\£xiAi]^ Vi r(^)...rK) 



We can then group these propagators and apply 
twice a Mellin-Barnes transformation 
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The Mellin-Barnes integral is defined with a con- 
tour parallel to the imaginary axis, which is 
placed at a value 7, such that it separates the 
poles of r(— a) from the poles of T(v + a). Since 
we have a sum of three propagators, we have to 
apply the Mellin-Barnes integral twice and ob- 
tain: 
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The term r(e, <r, r) is an abbreviation for a 
fraction of T-functions of the arguments E,Ui,a 
and t. The three-point function depends on the 
momenta of the two-point function via the terms 
with exponents a and r only. Combining this 
result with the result of the two-point function, 
which simply is a fraction of two T-functions, 
shifts their arguments by — r, —a respectively: 

T(u u -D/2) T(-a-T + v u -D/2) 



T(vu) 



T(-a 



In this way, we obtain for the two-loop two-point 

, ,. f (2,5) 

function 1 : 



f (2,5) 



71+zoo 

da 



72+ioo 
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where T'(e, z^, a, r) now consists of T-functions 
stemming from the result of the two-point and 
the three-point function. c(T) is a factor, which 



Table 1 

The first four Mellin moments for graphs I to VII, using M. Czakon's MB package. All v. L — 1 except for 



lb: v 1 = 2. 
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II 


-0.24999 £~ 
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-0.20833 £~ 
+0.46967 
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contains constants and other T-functions not de- 
pending on the arguments a and r. Closing the 
contour at infinity and collecting the residues of 
the T-functions inside the integration area via: 



we are left 



(-l) r 

res |r(— x + a), x — a + n] = — 

n! 

with a result which consists in general of (infinite) 
sums of T-functions: 



*(2,5) . . v-^ v-^ r(n ± h)T(n ± b 2 ) 

1 = nlT (n±b 3 ) 

r(j± Cl )r(j±c 2 ) 
j\ r(j ± c 3 ) 

r(n + j±a 1 )Y{n + j±a 2 ) 
T(n + j ± a3)T(n + j ± 04) 

These sums can be calculated with the computer 
libraries nestedsums [ 11 (C++) or XSummer [ I12j 
(FORM) . 

3. Operator Matrix Elements 

We now turn to the calculation of the diagrams 
with operator insertion. The previous calculation 
is changed due to the emergence of an additional 
numerator structure. As an example, let us con- 



sider graph I of Figure ^ 

The two-point function now carries an operator 
insertion, while the three-point function is unal- 
tered. To check for the change of the two-point 
function by the insertion, we first apply again the 
Feynman parameterization to it: 
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Throughout the calculation, we have to shift 
k\ — > k\ + x 2 p, leading to a corresponding shift 
in the numerator (A.ki) N — ► (A.(k\ + x 2 p)) N , 
which can be expressed via a binomial sum to be: 

N 

(A.k 1 +x 2 A.p) N = ( A (&-ki) l (x 2 A.p) N - 1 
1=0 ^ ' 

Since A is a light-cone vector and hence A 2 = 
0, we find that all integrals with (A.fci)' vanish, 



Table 2 

The first four Mellin moments for graphs I to VII. Wj = 1; lb: v\ = 2. 
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except the one for I = 0. This leads to a result 
for the two-point function, which is very similar 
to the original one, but which now contains the 
Mellin parameter N: 
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Here, we use = v\ + etc. Inserting the 
three-point function, we obtain for this graph I: 
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Analogously, we built the Mellin-Barnes inte- 
grals for the remaining six graphs and used the 
mathematica package MB by M. Czakon [E3j to 
numerically produce the results for the first few 
Mellin moments, given in Table 1. They serve as 
a check for our analytic result. 

To derive the analytic results, we continued 
from here and used relations like 

1 f + °° 

/ dsT(a + s)T(b + s)T(d - a - b - s) 

r(e - c + s)r(-s) 

x r(e + s) 

_ r(e - c)T(a)T(b)T(d - a)T(d - b) 
~ ne)T{d) 

x 3 F 2 [a, b, c; d, e; 1] , 

and other relations which are more general than 
Barnes formulae [ I13j . For the second integration 
we usually applied the Residue Theorem and in 
this way obtained double sums that contain the 
symbolic parameter N: 
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Table 3 

The analytic results for graphs I to VII for general values of N, with all vi = 1, lb: v\ — 2. 
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x [r(e)r(l-e)x 

r(j + 1 - 2e)r(j + 1 + £ ) T(fc + j + 1 + TV) 
r(j + 1 - e)T(j + 2 + N) T(k + j + 2) 
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Sums like these containing a symbolic parameter 
TV cannot be done via nestedsums or Xsummer, 
except for the most simple cases, fixing TV. For 
other more complicated cases, however, we had 
to use special transformations. As a computer 
algebra system we used MAPLE. For the first four 
Mellin moments the results are shown in Table 2 
and agree with the numeric results obtained by 
the use of MB. The graphs IV and V are, for all val- 
ues of Vi = 1, related to the graph Ia, respectively 
II, by a simple factor [1 + (—1)^] due to the fact 
that the operator insertion is located on a three- 
vertex with an on-shcll external line with p 2 = 0. 



Diagrams VI and VII require special treatment. 
The corresponding integrals are evaluated most 
effectively first transforming them into general- 
ized hypergeometric functions [ I13j . For fixed val- 
ues of TV again these analytic results agree with 
the numerical values presented in Table 1. 

Finally, we derive the complete analytic result 
for general values of TV. To obtain this, we had to 
extensively use algebraic and analytic relations 
to convert the sums obtained in evaluating the 
Mellin-Barnes integrals, into expressions contain- 
ing (nested) harmonic sums and related objects in 
intermediary steps. The final results depend on 
harmonic sums only and are summarized in Ta- 
ble 3. From this form the analytic continuation 
to complex values of TV [171 1141 ITS] needed in data 
analyzes 1 can be performed directly. Although 
the diagrams I- VII are the most sophisticated in 
the 2-loop problem under consideration, the re- 



1 A fast numerical analytic continuation of the heavy flavor 
Wilson coefficients [ 4 was given in f llbl . 



suits turn out to be very simple and much of the 
complexity of the Wilson coefficient stems from 
other, more simple structures. 

In summary, we have shown that by use of 
Mellin-Barnes integrals and direct representa- 
tions through generalized hypergeometric func- 
tions, massive two-loop five-propagator integrals 
containing operator insertions can be calculated 
in a very effective way in terms of nested har- 
monic sums. The corresponding expressions have 
a rather simple structure and can be expressed 
even in terms of only single harmonic sums. The 
intcgration-by-part method, on the other hand, 
allows to treat diagrams of lower complexity (4- 
propagator integrals), but leads to results with a 
much more complex structure. 
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